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Method jeg -For processing processing Rr aw Mmagnetic 
Resonance Ddata, Mmagnetic Rresonance Tfeomography 
Uttnit Aand Mraagnotic R*esonance Sepectroocopy Utmit 

[0001] The present application hereby claims priority 
under 35 U.S.C. §119 on German patent application 
number PS 10247150,9 filed September 10, 2002, the 
entire contents of which are hereby incorporated herein 
by reference. 

Field of the Invention 

f00021 The invention generally relates to a method for 
postprocessing raw magnetic resonance data. The 
invention also generally relates to a magnetic 
resonance tomography unit or a magnetic resonance 
spectroscopy unit which is respectively matched to such 
a method for postprocessing raw magnetic resonance 
data . 

Background of the Invention 

f00031 The aim of processing raw magnetic resonance data, 
which can originate both from magnetic resonance 
spectroscopy units and from magnetic resonance 
tomography units, is to extract the medically relevant 
information from the raw magnetic resonance data ais 
well as possible. In the case of magnetic resonance 
tomography, the aim is accordingly a picture with high 

resolution of detail and low noise^ T- — fehat is to say^_ 

the aim is for a picture with high edge sharpness and a 
high signal to noise ratio (SNR) . 
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f 00041 In the case of magnetic resonance- spectroscopy, 
the aim is to obtain from the raw data a linear 
spectrum with high resolution and with a high SNR. The 
quality in the processing of raw magnetic resonance 
data is thus firstly in the edge sharpness and image 
sharpness, i.e. in sharp pictures with high contrast^ 
— and a Secondly , it is in the production of a high SNR, 
so that the essence of the picture or of the spectrum 
is highlighted. 

fOOOSI In the case of the various approaches to achieving 
this aim, . it is usually not possible to increase the 
image sharpness and the SNR simultaneously. This 
drawback has a direct effect on the quality of the 
magnetic resonance tomography and spectroscopy. To 
improve the resolution of detail, radio- frequency 
"blue" noise is superimposed on a picture, for example. 
This produces a subjective increase in the resolution 
of detail. A likewise subjective improvement in the 
resolution is brought about in a similar manner by the 
application of high pass filters to the image data. 
This also amplifies the image noise, however. 

[00061 To improve the SNR in the case of magnetic 
resonance pictures, known filters (such as the Hanning 
filter, the Fermi filter or the cosine filter in the 
frequency domain) are applied to the raw magnetic 
resonance data. The action of such filters is known 
from the literature, e.g.: F.J. Harris. Proc . IEEE, 
Vol. 66, No 1 (1978): "On the Use of Windows for 
Harmonic Analysis with the Discrete Fourier Transform" , 
It is also possible to apply postprocessing filters in 
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the space domain, for example a mean filter, median 
filter or ARMA filter. 

[00071 In addition, US patent 4,463,375 discloses a 
method for medical image processing. This method 
specifies a reduced-noise version of a first processed 
image which has been obtained from a multiplicity of 
measurements in a multipictiire system, e.g. using a 
computer tomography unit. The method first involves 
producing a second image with a high SNR from the 
multiplicity of measurements . The first image is then 
processed using a filter which reduces the noise. The 
second image is processed using a filter which is 
complementary to this filter. The weighted combination 
of the two images produces the reduced-noise version of 
the first processed image. 

SUMMARY OF THE INVENTION 

[0008] An embodiment of t^ he invention « — booed — — fe^e 
includes an object of specifying a method for 
processing raw magnetic resonance data and a magnetic 
resonance tomography unit and a magnetic resonance 
spectroscopy unit which produce pictures or spectra 
with high edge sharpness without impairing the image 
noise, i.e. the SNR, in the process. 

r00091 For the method, an embodiment of the invention 
achieves thio — an object by virtue of a method having 
the method features below being carried out in order to 
postprocess raw magnetic resonance data. In a first 
step, the raw magnetic resonance data are filtered 

using a first filter-; T^-hese filtered data are 
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Fourier transformed in a second step . T- — the absolute 
value for the Fourier transformed data is then formed 
in a third step and in this way a first magnetic 
resonance signal is obtained. In a fourth step, the 
unfiltered, original raw magnetic resonance data are 
Fourier transformed . T- — and t his is followed in the 
fifth step by second absolute value formation for the 
unfiltered Fourier transformed data. Finally, weighted 
combination of the two magnetic resonance signals is 
performed and in this way the postprocessed magnetic 
resonance signal is obtained. 

^ — other — wordo , £The method comprioo - o — includes the 

weighted combination of two magnetic resonance signals, 
the magnetic resonance signals having been obtained 
from the same raw magnetic resonance data once with and 
once without filtering. 

[00101 * P Th rough suitable choice of the filter and of 
the parameters for the weighted combination, a picture 
which is optimized in terms of the representation 
wanted from a medical point of view is extracted. 

rOOin One advantage of the method is that the SNR is 
significantly improved following the weighted 
combination of the two magnetic resonance signals and 
that an edge - emphasized magnetic resonance signal is 
produced at the same time. In comparison with the 
method from US patent 4,463,375, the involvement for 
this method is significantly reduced . This occurs— 
since it is firstly based only on the raw data from a 
single measurement, and since it is secondly 
independent of an image with a high SNR which is 
produced from a plurality of measurements. 
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f00121 In one particularly simple extension of the 
method, the raw magnetic resonance data are filtered 
using a second filter before the second magnetic 
resonance signal is produced. This allows the second 
magnetic resonance signal to be influenced and thus 
allows it to improve the weighted combination. The use 
of the unfiltered raw magnetic resonance data in the 
original method corresponds to the application of the 
identity filter and is faster to execute. 

r0013]_ A weighted combination within the context of an 
embodiment of the invention is in this case a 
mathematical combination of the absolute values of the 
two magnetic resonance signals. These have the same 
local coding and represent the echo signal from a 
magnetic resonance measurement, for example in a 
one-dimensional or multidimensional pixel structure, 
following the Fourier transformation. In this case, the 
combination, which is usually addition, is spatially 
consistent, i.e. the same spatial areas are combined. 

[00141 By way of example, combination takes place pixel 
by pixel. The weighting is effected using a weighting 
factor. For a spatial area of one of the two magnetic 
resonance signals, this weighting factor depends on the 
information content, e.g. the intensity, of the same 
spatial area in the other magnetic resonance signal. By 
way of example, the weighting factor can have a linear 
dependency on the information content or can have any 
other beneficial nonlinear dependency on the 
information content of the spatial area. The weighting 
factor can be produced, by way of example, by a 
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mathematical step function which sets the weighting 
factor to a small value or zero below a previously 
stipulated threshold value for the absolute value of 
the magnetic resonance signal in question. 

rQ015L One advantage of weighted combination is that the 
SNR is significantly improved following the weighted 
combination of the two filters . Further, and tha ^ this 
data weighting produces an edge - emphasized profile in 
the magnetic resonance signal. 

[00161 The method relates to the postprocessing of raw 
magnetic resonance data from a single measurement and 
is not dependent on whether the raw data belong to a 
single picture or whether they belong to a picture from 
a series of pictures which show a development over 
time, for example. The raw magnetic resonance data can 
be raw data from a measurement in a magnetic resonance 
spectroscopy unit or else raw data from a picture in a 
magnetic resonance tomography unit. 

[00171 Accordingly, an embodiment of . the invention 

achieves a4rhe n object for a magnetic resonance 

tomography unit by moono — use of a magnetic resonance 
tomography unit which is matched to a method for 
postprocessing raw magnetic resonance data with the 
method features above. In this case, the raw magnetic 
resonance data are referred to as a spin echo signal or 
an echo signal, for example. The magnetic resonance 
signals correspond to the one-dimensional or 
multidimensional magnetic resonance pictures, for 
example . 
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f00l81 An embodiment of t ghe invention achieves an-&i*e 
object for a magnetic resonance spectroscopy unit by 
mcano use of a magnetic resonance spectroscopy unit 
which is matched to a method for postprocessing raw 
magnetic resonance data with the method features above. 
In this case, spectroscopic data, e.g. the FID (Free 
Induction Decay) signal, are postprocessed . 

r00191 In one particularly beneficial embodiment of the 
method, the filters are applied to the raw magnetic 
resonance data after demodulation thereof. This allows 
the filters to be applied in their customary form of 
representation in the frequency domain. Applying the 
filters after demodulation has the advantage that the 
filters can be very narrowband. 

r00201 In one particularly advantageous embodiment, the 
first filter is a low pass filter which filters rapid 
changes out of the raw data and thus suppresses radio- 
frequency oscillations. By way of example, the low pass 
filter can be a type of Hanning filter. 

r00211 In another embodiment the second filter is a high 
pass filter which passes only the rapid changes in the 
signal, which are also caused by noise, for example. 
Magnetic resonance signals which have been high pass 
filtered show an excessive increase in rapidly rising 
edges . 

[00221 In one particularly advantageous embodiment which 
uses self-weighting for the combination, the weighted 
combination is effected as follows: 



7 



OCT-0S-2003 10:42 HDP 7033903020 P. 10 



New Fatent Application 
Docket No. 32860-0006I9/US 



C = A + /1 



B 



[Equation l] 



r00231 In this case, A and B are magnetic resonance 
signals which are subjected to weighted combination by 
mcano of -low pass filtering (A) and high pass filtering 
(B) , for example. A^ ax is the maximum of the magnetic 
resonance signal A. The two parameters X and k 
determine the contribution level of the magnetic 
resonance signal B to the postprocessed magnetic 
resonance signal C. The quotient of B and A^ brings 
about amplified correction of the pixels which have a 
high absolute value in the magnetic resonance signal S. 
The values of X and k need to be ascertained 
empirically on the basis of their influence and their 
optimum effect on the magnetic resonance signal. In 
this case, the values of X and K are preferably between 
1 and 3 . 

f00241 One advantage of the quotient which normalizes B 
to A max is that an uncontrolled signal contribution for 
B is avoided in regions in which there is actually 
nothing to correct or to improve. 

[00251 T he weighted combination of suitably filtered 
magnetic resonance signals has the advantage that, in a 
single magnetic resonance measurement, the noise in the 
resultant magnetic resonance signal is lower than with 
linear addition or subtraction of the two filtered 
magnetic resonance signals. 
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r00261 In particular embodiments, the method from the 
invention can be applied to raw magnetic resonance data 
which cover a one-dimensional or multidimensional space 
which is to be examined. Depending on the application, 
one-dimensional sections, two-dimensional sectional 
views or three-dimensional volume representations are 
shown, for example, in magnetic resonance technology. 
The raw magnetic resonance data contain the dimension 
of the respective picture. 

[00271 In one advantageous embodiment, the dimensionality 
of the Fourier transformation performed during the 
method is chosen such that the dimension of the raw 
magnetic resonance data is retained. Th± a mcano that us^ 
the method can be applied, by way of example, to two- 
dimensional magnetic resonance sectional views or to 
three-dimensional volume representations. Even in the 
case of magnetic resonance spectra, one-dimensional or 
multidimensional representation, for example in three 
dimensions, is possible. 

[00281 In one particularly advantageous embodiment of the 
method, two locally coded magnetic resonance signals 
are combined such that the contribution of one magnetic 
resonance signal to the weighted combination is formed 
by multiplying this magnetic resonance signal by a 

weighting factor-; w^rfeh fe. T he weighting factor 

depends**^ on the other magnetic resonance signal such 
that it is greater for a large magnetic resonance 
signal than for a small magnetic resonance signal. 

f 00291 In one specific embodiment^ the weighting factor 
has a nonlinear dependency on the absolute value of a 
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magnetic resonance signal. At least in one dimension of 
the raw magnetic resonance data, the same nonlinear 
dependency is used in order to have the same effect on 
the image sharpness, for example, over all the pixels 
in this dimension. 

[00301 In one particular embodiment of the method, the 
raw magnetic resonance data are processed using a 
plurality of filters before the differently filtered 
magnetic resonance signals are subsequently combined by 
means of weighted combination to form a single 
postprocessed magnetic resonance signal. 

BRIEF DESCRIPTION OF THE DRAWINGS 

f0031] A plurality of exemplary embodimen ts of the 
method, of the magnetic resonance tomography unit and 
of the magnetic resonance sp ectroscopy unit based on 
the invention will become more fully understood from 
the detailed description of preferred embodiments given 
hereinbelow and the accompanying drawings, which are 
given bv way of illustration onl y and thus are not 
limitative of the present invention, and wherein: 

A plurality e~§ e xemp la ry embodiment o &4 fe^e 



method, — of the magnetic resonance tomogr a phy unit D T> ri 

e£ the magnetic reoonancc — s pectroscopy unit baaed on 

the invention arc explained in more detail below wren 
reference to-figurco 1 to 8 , — in which ;- 

figure 1 shows a flowchart of an exemplary embodiment 
of the method which has been implemented in a 
magnetic resonance tomography unit, 
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figure 2 shows a flowchart of an exemplary embodiment 
of the method in which two filtered magnetic 
resonance signals are subjected to weighted 
combination , 

figure 3 shows a numerical simulation of the 
processing of raw magnetic resonance data 
using a spin echo signal which is brought 
about by a square object in the image 
gradient of a magnetic resonance tomography 
unit, 

figure 4 shows a numerical simulation of the 
processing when a low pass filter is applied 
to the spin echo signal from figure 3, 

figure 5 shows a numerical simulation of the 
processing when a high pass filter is applied 
to the spin echo signal from figure 3 , 

figure 6 shows an abscissa detail from figures 3 to 5 
to illustrate the action of the method using 
the edge- emphasized profile which is achieved 
by the weighted combination, 

figure 7 shows an enlargement of figure 6 in the area 
of an edge of the square object, 

figure 8 shows a numerical simulation of the 
processing of raw magnetic resonance data for 
the case of magnetic resonance spectroscopy 
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using an FID signal which has been formed by 
the sum of three Lorentz lines, and 

figure 9 shows a comparison of the magnetic resonance 
spectrum obtained by weighted combination 
with the low pass filtered magnetic resonance 
spectrum, both magnetic resonance spectra 
having been normalized to the same noise. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[00321 Figure 1 schematically shows the flow of the 
method and the necessary components which are required 
within a magnetic resonance tomography unit 1 for the 
flow of the method. In a magnetic resonance tomography 
unit 1 having a conventional magnet and gradient system 
(not shown explicitly) used for local coding, a radio- 
frequency transmitter 3 is used to irradiate an 
examined object 5 with radio- frequency pulses. The 
electromagnetic signals emitted by the examined object 
5 are received using a radio-frequency receiver 7. The 
output signals from the radio-frequency receiver 7 form 
the raw magnetic resonance data in the form of a spin 
echo signal 9 . 

[00331 These data are postprocessed within a data 
processing installation 11. The latter subjects them to 
weighted combination with one another either in 
unfiltered form or in a form filtered using various 
filters 13 in a weighting unit 15. The result is a 
postprocessed magnetic resonance signal 17 in the form 
of a magnetic resonance picture. 
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r00341 Figure 2 shows a flow chart of an exemplary 
embodiment of the method in which two filtered magnetic 
resonance signals are subjected to weighted combination 
in the data processing installation 11. In this case, 
the spin echo signal 9 is first demodulated in a 
demodulation unit 19, so that it is available in the k 
domain with the dimension of the sampling. Next, the 
spin echo signal 9 is supplied to two data processing 
units 21 having a respective filter. There, the spin 
echo signal is filtered once using the low pass filter 
processor 23 and the other time using the high pass 
filter processor 25. 

fQ0351 Following the subsequent Fourier transformations 
27 and the absolute value formations 29, one firstly 
obtains the low pass filtered magnetic resonance 
picture 31 and the high pass filtered magnetic 
resonance picture 33, respectively. These two magnetic 
resonance pictures are combined in the weighting unit 
15 by weighted addition to form a single postprocessed 
magnetic resonance picture 35, The weighted addition is 
carried out by inputting parameters 37 using an input 
unit 39. The parameters 37 determine the weighting of 
the pixel 41 in the high pass filtered magnetic 
resonance picture 33 on the basis of the same pixel 43 
in the low pass filtered magnetic resonance picture 31. 
With a suitable choice of weighted addition and its 
parameters, the postprocessed magnetic resonance 
picture 3 5 has both a high SNR and great edge 
sharpness . 

[00361 W hen the method is used with a magnetic resonance 
spectroscopy unit, the flowchart is comparable with 
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that in figures 1 and 2 . The only difference is that 
raw data from magnetic resonance spectra are filtered 
and subjected to weighted addition by means of the 
method in order to achieve a higher resolution and a 
greater SNR. 

[00371 Figures 3 to 5 show the results of a numerical 
simulation of a spin echo signal for a square object in 
the field gradient of the magnetic resonance tomography 
unit 1. The simulation relates to a section through a 
square object and illustrates the effect of the method 
on the object's edge sharpness. The object is in the 
area of the 128 pixels in the simulation. 

[00381 The figures respectively show the real and 
imaginary parts 51 and 53 of the echo signal after a 
filter has been applied to the echo signal. They show 
the filter function and the complex-value spin echo 
signal above a k s axis in the k domain. Various filter 
functions clarify the action of the filters on the echo 
signal in the K domain. In addition, figures 3 to 7 
show the result of Fourier transforming the filtered 
complex-value echo signal with subsequent absolute 
value formation. In this way, a simulated measurement 
result is obtained from a magnetic resonance tomography 
unit, in this case a one-dimensional section through 
the square object. The sectional view of the 
measurement result is plotted above the corresponding 
spatial axis. 

f00391 Figure 3 shows the real part 51 and the imaginary 
part 53 of the unfiltered echo signal and also the 
result of the Fourier transformation and of the 
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absolute value formation for the unfiltered echo 
signal, in this case referred to as identity simulation 
result 57. In the simulation, no filter 55 was applied 
to the echo signal in this case. The oscillations in 
the identity simulation result 57 in the area of the 
rising and falling edges can be seen clearly. These 
oscillations limit the edge sharpness in the simulated 
(unpostprocessed) measurement- 

r00401 Figure 4 shows the influence of the low pass raw 
data filter 59 on the spin echo signal. The low pass 
raw data filter 59 has been applied in the form of a 
Hanning filter F* to the spin echo signal, the Hanning 
filter F H having the following form in the 
one-dimensional k domain; 

F TM =F H =l/2x[l+co S (27lO-A72)/N)] [Equation 2] 

fOQ411 In this case, N is the number of sampling points 
in the picture in one dimension. The filter modulation 
follows a cosine function whose maximum is shifted by 
N/2, that is to say to the center of the k domain. The 
index j describes the jth sampling point. The filter 
suppresses the marginal area of the k domain with lower 
intensity in the echo signal. This can clearly be seen 
from the reduced oscillation amplitude of the filtered 
spin echo signal in the marginal area. The low pass 
simulation result 61 in figure 4 shows a rounded edge 
for the square object. This can also be seen in figures 
5 and 6 . 
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fQQ421 in figure 5, the simulation is carried out for a 
high pass raw data filter 63, with the high pass filter 
in the k domain being a modulated Hamming filter F H 
provided with an offset: 

Fhh = ( Fh +0j ^)xcos n M;-N/2)/N) [Equation 3] 

f 0 Q431 in this case, the modulation of the Hanning filter 
F„ takes place with a cosine function to the power of n, 
with the minima of the cosine function being at the two 
ends and in the center of the k domain. The additional 
minima at the ends have the advantage that they cause 
no additional artefacts in the Fourier transformation. 
Thus, by way of example, the "truncation artefact" in 
the Fourier transformation is suppressed. 

roQ441 in this case, the offset parameter of the 
modulated Hanning filter has the value 0.2 and the 
power n of the cosine modulation has the value 1. The 
high pass raw data filter 63 amplifies the wings of the 
echo signal in the k domain. In the real part 51 and in 
the imaginary part 53 of the filtered echo signal, it 
is possible to see the amplified amplitude of the echo 
signal oscillations in the marginal areas and a 
suppression in the central area. The high pass 
simulation result 65 shows a marked overshoot in the 
area of the edges of the square object. 

fOQ45l Figure 6 shows a detail from the simulation 
results in the area of the pixels 35 to 55, i.e. in the 
area of the rising edge. The oscillations in the 
identity simulation result 57 can clearly be seen. The 
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edge of the square object is rounded in the case of the 
low pass simulation result 61. 

rOQ461 The high pass simulation result 65 shows a 
pronounced signal overshoot right in the area of the 
rise. Using the weighted addition of the two filtered 
simulation results, it is now possible to achieve an 
edge - emphasized postprocessing result 67. To this 
end, a weighted combination has been performed in line 
with the equation 1 already indicated above. In this 
case, the low pass simulation result 61 corresponds to 
the magnetic resonance signal A and the high pass 
simulation result 65 corresponds to the magnetic 
resonance signal B . 



r0 0471 The parameters X and K are set such that a maximum 
sharpness is produced in the edge area. In this case, X 
and K are usually between the values 1 and 3 . The level 
obtained for the weighting of the contribution of B 
depends on the ratio of the absolute value in the 
respective pixel in B to the maximum value A ma x of the 
pixels in A, i.e. the higher the value of the high pass 
simulation result 65 in a pixel, the higher its 
contribution. X then performs a type of linear 
weighting, while k weights the influence of this 
quotient as a power of the ratio of B and A, that is to 
say indicates how greatly the correction of A is 
determined by the ratio of magnitudes for B with 
respect to A. The parameters X and k have been set to 
Che values 0.6 and 2.3, respectively, in figures 6 and 
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r00481_ Besides the edge sharpness, the SNR of the 
postprocessing result 67 is also increased, and is now 
comparable to the SNR of the low pass simulation result 
61. This can be seen from the suppression of the signal 
before the rising edge, i.e. there is a low base 
contribution outside the object. 

f 00491 Figure 7 shows a further enlarged detail from the 
simulated measurement results. The representation 
focuses on the edge area of the square object in the 
rea of pixels 42 to 55 with signal levels between 10 
and 13. It can clearly be seen how the greatly 
oscillating identity simulation result 57 is freed from 
the radio-frequency oscillations using the low pass 
filter. Next, the weighted addition fills the comer 
area of the square object on account of the high signal 
absolute values in the high pass simulation result 65. 
The oscillations in the identity simulation result 57 
are also called "ringing". This artefact, which is also 
called truncation artefact, is suppressed using the 
method of the invention. 
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r00501 Figure 3 shows the flow of the method applied to 
raw magnetic resonance data from a simulated magnetic 
resonance spectroscopy measurement in which an FID 
signal 71 is measured which is formed by the sum of 
three Lorentz lines (73, 75, 77). Figure 3 shows the 
real parts of the FFT spectra involved in the method. 
In the simulation, a noise contribution has been added 
to the synthesized FID signal 71. Low pass filtering of 
the FID signal 71 results in the low pass filtered 
spectrum 79, and high pass filtering results in the 

IS 



□CT-06-2003 10:44 



HDP 



703390.3020 P. 21 



« ' New Patent Application 

Docket No. 32860-000619/US 

high pass filtered spectrum 31. Following the weighted 
combination of the two filtered spectra, the 
postprocessed magnetic resonance spectrum 83 is 
obtained. It is possible to see markedly reduced noise 
and high resolution of the three Lorentz lines on 
account of the use of the method of the invention. 

[00511 In figure 9, the postprocessed magnetic resonance 
spectrum 83 is compared with the low pass filtered 
spectrum 79. To clarify the influence of the method, 
the two spectra have been normalized to the same noise 
in this case. The narrower line widths and a 
considerable increase in the contrast in the 
postprocessed magnetic resonance spectrum 83 are 
clearly evident. 

f00521 Exemplary embodiments being thus described, it 

will be obvious that the same may b e varied in many 
ways. Such variations are not to be regarded as a 
departure from the spirit and scope of the present 
invention. and all such mo difications as would be 
obvious to one skilled in th e art are intended to be 
included within the scope of the f ollowing claims. 
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